Fractional variational problems 
depending on indefinite integrals* 

O | Ricardo Almeida Shakoor Pooseh Delfim F. M. Torres 

ricardo.almeida@ua.pt spooseh@ua.pt delfim@ua.pt 



O 



o 



>< 



Department of Mathematics, University of Aveiro, 3810-193 Aveiro, Portugal 



Abstract 



^) [ We obtain necessary optimality conditions for variational problems with a Lagrangian 



depending on a Caputo fractional derivative, a fractional and an indefinite integral. Main 
results give fractional Euler-Lagrange type equations and natural boundary conditions, which 
i-SZj . provide a generalization of previous results found in the literature. Isoperimetric problems, 

problems with holonomic constraints and depending on higher-order Caputo derivatives, as 
well as fractional Lagrange problems, are considered. 
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1 Introduction 

CM . In the 18th century, Euler considered the problem of optimizing functionals depending not only 



on some unknown function y and some derivative of y, but also on an antiderivative of y (see |19j). 
Similar problems have been recently investigated in |24j . where Lagrangians containing higher- 
order derivatives and optimal control problems are considered. More generally, it has been shown 
that the results of [24] hold on an arbitrary time scale [30]. Here we study such problems within 
the framework of fractional calculus. 



Roughly speaking, a fractional calculus defines integrals and derivatives of non-integer order. 



Let a > be a real number and n € N be such that n — 1 < a < n. Here we follow [5] and [2l)ll3T] . 
Let / : [a, b] — > R be piecewise continuous on (a, b) and integrable on [a, b]. The left and right 
Riemann-Liouville fractional integrals of / of order a are defined respectively by 



1 r x 1 rb 



=pfr/ {x-tT-'f(t)dt and a I«f{x) = -j-l (t - x)"' 1 f(t)dt. 
J a "(a) J x 

Here T is the well-known Gamma function. Then the left a D" and right X D^ Riemann-Liouville 
fractional derivatives of / of order a are defined (if they exist) as 



T(n - a) dx n J a 

and 
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The fractional derivatives (JXJ) and @ have one disadvantage when modeling real world phenomena: 
the fractional derivative of a constant is not zero. To eliminate this problem, one often considers 
fractional derivatives in the sense of Caputo. Let / belong to the space AC n ([a, b];M) of absolutely 
continuous functions. The left and right Caputo fractional derivatives of / of order a are defined 
respectively by 

1 



and 



1 ' b 



1 D%f{x) = - / {-l) n {t-x) n - a - l f n \t)dt, 

1 [n — a.) 



These fractional integrals and derivatives define a rich calculus. For details sec the books 26 31.39 . 
Here we just recall a useful property for our purposes: integration by parts. For fractional integrals, 

r b pb 

9(x) ■ J£f(x)dx = / f(x) ■ x I£g(x)dx 

J a 

(see, e.g., [351 Lemma 2.7]), and for Caputo fractional derivatives 

r b r b n-i b 

/ g(x) ■ C a D a x f{x)dx = / f(x) ■ x DZg(x)dx + £ \ x D a b + ^ n g{x) ■ f {n ~^\x) 

J a J a j=0 

(see, e.g., [3l Eq. (16)]). In particular, for a S (0, 1) one has 

b g(x) ■ C a D a J(x)dx = f f{x) ■ x D«g(x)dx + [j^gix) ■ f{x)] \ . (3) 

J a 

When a — > 1, ^5" = 4-, x D b = — x ll~ a is the identity operator, and (J3]) gives the classical 
formula of integration by parts. 

The fractional calculus of variations concerns finding extremizers for variational functionals 
depending on fractional derivatives instead of integer ones. The theory started in 1996 with the 
work of Riewe, in order to better describe non-conservative systems in mechanics |37il38) . The 
subject is now under strong development due to its many applications in physics and engineering, 
providing more accurate models of physical phenomena (see, e.g., [Hl9l lT2l[Lll[T? , 18,20 ,21,33,55]). 
With respect to results on fractional variational calculus via Caputo operators, we refer the reader 
to HSIIimCMiniGIl] and references therein. 

Our main contribution is an extension of the results presented in [3][53] by considering La- 
grangians containing an antiderivative, that in turn depend on the unknown function, a fractional 
integral, and a Caputo fractional derivative (Section [5]) . Transversality conditions are studied in 
Section |31 where the variational functional J depends also on the terminal time T, J(y,T), and 
where we obtain conditions for a pair (y, T) to be an optimal solution to the problem. In Sectionf?] 
we consider isoperimetric problems with integral constraints of the same type as the cost func- 
tionals considered in Section [2] Fractional problems with holonomic constraints are considered in 
Section [5j The situation when the Lagrangian depends on higher order Caputo derivatives, i.e., it 
depends on ^D" k y(x) for ak £ (k — 1, fc), fc G {1, . . . is studied in Section|5J while Section[7] 
considers fractional Lagrange problems and the Hamiltonian approach. In Section |S] we obtain 
sufficient conditions of optimization under suitable convexity assumptions on the Lagrangian. We 
end with Section [9j discussing a numerical scheme for solving the proposed fractional variational 
problems. The idea is to approximate fractional problems by classical ones. Numerical results for 
two illustrative examples are described in detail. 



2 The fundamental problem 

Let a € (0, 1) and (3 > 0. The problem that we address is stated in the following way. Minimize 
the cost functional 

J(y)= [ L(x,y(x)CD«y(x), a I%y(x),z(x))dx, (4) 
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where the variable z is defined by 

z(x)= f l(t,y(t)CD?y(t), a I?y(t))dt, 

subject to the boundary conditions 

y(a) = y a and y(b) = y b . (5) 

We assume that the functions (x, y, v, w, z) L(x, y, v, w, z) and (x, y, v, w) — > l(x, y, v, w) are of 
class C , and the trajectories y : [a, b] — > K are absolute continuous functions, y £ AC([a, b];M), 
such that ^D^y(x) and a Igy(x) exist and are continuous on [a, b]. We denote such class of 
functions by FQa, 6];M). Also, to simplify, by [•] and {•} we denote the operators 

[y](x) = {x,y{x), c a D«y{x), a lPy{xlz(x)) and {y}(x) = (x, y(x), c a D«y{x), a I*y{x)). 

Theorem 1. Let y £ J-([a,b];M) be a minimizer of J as in subject to the boundary conditions 
([5]). Then, for all x £ [as, b], y is a solution of the fractional equation 

f y [y\{ X ) + ^ + m ig (^ MW ) + [f z [y\m- d ^{vm 



dL 

dz 



di 



[y](t)dt--{y}(x))+ x I^ 



dv 



dL 



01 



. )z [y](t)dt-—{y}(x)\=0. (6) 



Proof. Let h £ J-([a, b]\ K) be such that h(a) = = h(b), and e be a real number with |e| <C 1. If 
we define j as j(e) = J(y + eh), then j'(0) = 0. Differentiating j at e = 0, we get 

" V>>1 " -x)h(x) + ^[y}(x^D^h(x) + ^[y}(x) a I^h(x) 



a [dy 



dL 
"9z l 



dw 

J£{t/K*)M*) + ^{yW) c a D?h{t) + ^L{y}(t)jfh(i - j ,// 



dx = 0. 



The necessary condition © follows from the next relations and the fundamental lemma of the 
calculus of variations (cf., e.g., [?TJ p. 32]): 



^ ^[y]{x) c a D«h{x)dx = j x Dt{^[y]{x)\h{x)dx^ 

[y]{x)jth{x)dx = J x lg (f^[«) Hx)d. 



b dL 
dw 



dL 



y](x) 



x dl \ rb 

— {y}(t)h(t)dtj dx = 



dL 



dx d 



[y](t)dt 



x dl 



{y}(t)h(t)dt dx 



fa {[^ m)d l ^^ h ^ dx > 



a dy 



i or *«*) *-/.(-=/. &h or ° ? 



<9z 



a \ ./ re 



at 

9 2! 



(9/ 



^[y}(t)d&y}(x)) h(x)dx + 



dz 



dv 



x ± b 



" ^[y}(t)d&y}(x)) h(x) 



dz 



dv 



3 



and 



□ 



The fractional Euler-Lagrange equation © involves not only fractional integrals and fractional 
derivatives, but also indefinite integrals. Theorem Q] gives a necessary condition to determine the 
possible choices for extremizers. 

Definition 2. Solutions to the fractional Euler-Lagrange equation ([5]) are called extremals for J 
defined by Q. 

Example 3. Consider the functional 

J(y) = [ [(%D2v(x) - T(a + 2)xf + z(x)} dx, (7) 
Jo 

where a G (0,1) and 

defined on the set 
Let 
Then, 



*0=) = / (y(t)-t a+ ydt, 
Jo 

{y 6 7([0,1];R) : y(0) = and y(l) = 1} . 

y a (x)=x a+x , xe[0,l]. (8) 



%D2y a (x)=r(a + 2)x. 

Since J{y) > for all admissible functions y, and J(y a ) — 0, we have that y a is a minimizer of 
J. The Euler-Lagrange equation applied to gives 

m D?($DZy(x)-T(a + 2)x)+ [ ldt (y(x) - x a+1 ) = 0. (9) 



Obviously, y a is a solution of the fractional differential equation @. 

The extremizcr (JSJ of Example [3] is smooth on the closed interval [0, 1]. This is not always the 
case. As next example shows, minimizers of (|4"]|— ([5| are not necessarily C 1 functions. 

Example 4. Consider the following fractional variational problem: to minimize the functional 

J(V) = / \{oD a x y{x) - l) 2 + z(x)] dx (10) 

ye^([0,l];K) : y(0) = and y(l) = \ 

T(a + 1) 

where z is given by 

z{x) = fo{ m ^n^]) dL 

Since q Df:x a = T(a + 1), we deduce easily that function 



1 [a + 1) 
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is the global minimizer to the problem. Indeed, J(y) > for all y, and J(y) = 0. Let us see that 
y is an extremal for J . The fractional Euler-Lagrange equation ^ becomes 

2 x D?{°D«y{x) - 1) + jf 1 dt ■ 2 (y(x) - j^yy) = 0. (12) 



Obviously, y is a solution of equation {T 

Remark 5. The minimizer of Example^ is not differ entiable at 0, as < a < 1. However, 
17(0) = and $D%y~(x) = oD%y(x) = T(a + 1) /or any x € [0, 1]. 

Corollary 6 (cf. equation (9) of [2]). If y is a minimizer of 

J(y)= f L{x,y{x)° a D«y{x))dx, (13) 



.Sff 



6jec< to tte boundary conditions ([5j) 7 i/ien y is a solution of the fractional equation 



Proof. Follows from Theorem [T] with an L that does not depend on a I^y and z. □ 

We now derive the Euler-Lagrange equations for functionals containing several dependent 
variables, i.e., for functionals of type 

,6 

J(yi,..., y n ) = I L(x, yi{x), y n (x), ^D^y x (x), %D%y n (x), a I^y x (x), a l£y n (x),z(x))dx, 



(14) 

where n £ N and z is defined by 

z(x) = f l(t,y x (t), . . .,y n (t),%D?y x (t), . . .^Dfy n (t), a I^y x (t),.. . , a l£ y n (t))dt, 



subject to the boundary conditions 

yk(a)=y a ,k and y k (b) = y b ,k, k £ {l,...,n}. (15) 

To simplify, we consider y as the vector y = (yi, . . . , y n ). Consider a family of variations y + eh, 
where |e| <C 1 and h = (hi, . . . , h n ). The boundary conditions (IT51) imply that hk(a) = = hk(b), 
for k £ {!,... ,n}. The following theorem can be easily proved. 



Theorem 7. Let y be a minimizer of J as in (| 14[) , subject to the boundary conditions (|15p . Then, 
for all k £ {1, . . . , n} and for all x £ [a,b], y is a solution of the fractional Euler-Lagrange equation 

3 Natural boundary conditions 

In this section we consider a more general question. Not only the unknown function y is a variable 
in the problem, but also the terminal time T is an unknown. For T £ [a, b], consider the functional 

J(y,T)= f T L[y](x)dx, (16) 
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where 



[y](x) = (x,y(x)°D2y(x), a lgv(x),z(x)). 



The problem consists in finding a pair (y, T) € F([a, b]; R) X [a, b] for which the functional J attains 
a minimum value. First we give a remark that will be used later in the proof of Theorem [9] 

Remark 8. If (f) is a continuous function, then (cf. \31{ p. 46]) 

Mm x I*r a <t>(x) = 

for any a € (0, 1). 



Theorem 9. Let (y,T) be a minimizer of J as in (|16p. Then, for all x £ [a,T], (y,T) is a 
solution of the fractional equation 



-I 



(x) + X D% 



dL 
dv^ 



xDrp 



0)+^(^M(«))+ fTdL 

T ^[y](t)dt.^{y}( X y\ + x 4 



. h [y](t)dt-^-{y}(x) 



T dL , , . . , dl , , , . \ 



and satisfies the transversality conditions 



x ±r p 



dv 



<9z 



<9v 



T[y](T) = 0. 

Proof. Let ft- G -T([a, b]; R) be a variation, and let AT be a real number. Define the function 

i(e) = Jfo + e^T + eAT) 

with |e| <C 1. Differentiating j at e = 0, and using the same procedure as in Theorem [TJ we deduce 
that 



AT.L[y](T)- 



F)T ( BT 



z' 3 

T 



dL 



' " L [y](t)dt.^-{y}(x) 



dz 



dy 



x ±r p 



dL 
~dv~ [ 





T 




)h(x) 


+ 






a 





h(x)dx 
T 



" ' 5i -M(*)dt-#{v}(a:)U(«) 



02 



<9v 



0. 



The theorem follows from the arbitrariness of h and AT. □ 
Remark 10. If T is fixed, say T = b, then AT = and the transversality conditions reduce to 



l-a 



= 0. 



(17) 



Example 11. Consider the problem of minimizing the functional J as in (|10p. but without given 
boundary conditions. Besides equation (1121) . extremals must also satisfy 



[ x l{- a (^D^y(x)-l)] o = 0. 
Again, y given by (fTTj) is a solution of (|T2|) and (p~8|) . 



(18) 



G 



As a particular case, the following result of [2] is deduced. 

Corollary 12 (cf. equations (9) and (12) of [2]). If y is a minimizer of J as in (|13[) . then y is a 
solution of 

|blW+^(f[!/]W)=0 



and satisfies the transversality condition 

. ( 8L. . M 

= 0. 



Proof. The Lagrangian i in (1131) does not depend on a J^y and 2, and the result follows from 
Theorem [H □ 



Remark 13. Observe that the condition 

fr)T, \ 1 

= 

Jib 

is implicitly satisfied in Corollary (cf. Remark^. 



4 Fractional isoperimetric problems 

An isoperimetric problem deals with the question of optimizing a given functional under the 
presence of an integral constraint. This is a very old question, with its origins in the ancient 
Greece. They where interested in determining the shape of a closed curve with a fixed length and 
maximum area. This problem is known as Dido's problem, and is an example of an isoperimetric 
problem of the calculus of variations [41]. For recent advancements on the subject we refer the 
reader to [6H71H3E7] and references therein. In our case, within the fractional context, we state 
the isoperimetric problem in the following way. Determine the minimizers of a given functional 

,6 

J(y) = L(x,y(x)°BSy(x), a lj>y(x),z(x))dx (19) 



subject to the boundary conditions 

y(a) = y a and y(b) = y b (20) 
and the fractional integral constraint 



6 

Cr>«„ 



I(y)= G(x,y(x)^D«y(x), a lPy(x),z(x))dx= % 7 e K, (21) 

J a 

where z is defined by 

z{x)= [ l{t,y{t), c a D?y{t), a I?y{t))dt. 



As usual, we assume that all the functions (x, y, v, w, z) — > L(x, y, v, u>, z), (x, y, v, w) — > l(x, y, v, w), 
and (x, y, v, w, z) — > G(x, y,v,w, z) are of class C 1 . 

Theorem 14. Let y be a minimizer of J as in (|19|) . under the boundary conditions (I20p and 
isoperimetric constraint (|2ip . Suppose that y is not an extremal for I in J21|) . Then there exists 
a constant A such that y is a solution of the fractional equation 



where F = L — AG, for all x 6 [a, b] . 
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Proof. Let ei, e 2 € M. be two real numbers such that |ei| <C 1 and \e 2 \ <C 1, with t\ free and e 2 to 
be determined later, and let h\ and h 2 be two functions satisfying 



Define functions j and i by 
and 



h 1 (a) = h 1 (b) = h 2 (a) = h 2 (b) = 0. 
j(ei,e 2 ) = J{y + e\hi + e 2 h 2 ) 



i(ei, £2) = J(y + ei/11 + £2/12) - 7- 
Doing analogous calculations as in the proof of Theorem [TJ one has 



di_ 

de 2 



(0,0) J a 



dG r 



dy 



dG 



— [y](x)+ x D% — [y](x) + X I? — [y](x) 



dv 



dG r 



dw 



^[ y ](t)dt.^{y}(x) 



dG 



81 



()z [y](t)dt--{y}{x)\+ x I P b 



dG 



81 



h 2 (x) dx. 



By hypothesis, y is not an extremal for / and therefore there must exist a function h 2 for which 

?0. 



de 2 



(0,0) 



Since i(0, 0) = 0, by the implicit function theorem there exists a function e 2 (-), defined in some 
neighborhood of zero, such that 

i(ei,ca(ei))=0. (22) 



On the other hand, j attains a minimum value at (0,0) when subject to the constraint (|22p . 
Because Vz(0, 0) 7^ (0, 0), by the Lagrange multiplier rule |41[ p. 77] there exists a constant A such 
that 

V(j(0,0)-Ai(0,0)) = (0,0). 



So 



de 1 



- A 



(0,0) 



dt_ 



= 0. 



(0,0) 



Differentiating j and i at zero, and doing the same calculations as before, we get the desired 
result. □ 

Using the abnormal Lagrange multiplier rule |41| p. 82], the previous result can be generalized 
to include the case when the minimizer is an extremal of /. 

Theorem 15. Let y be a minimizer of J as in p9[) . subject to the constraints (|20p and (|2ip . 
Then there exist two constants Aq and X, not both zero, such that y is a solution of equation 



dK 



— [ y ]{ x ) +x D<g[—[y]{x 



dy 



OK, 



dv 



dK, ,, A [ b dK, „ . , 81 ... . 



b dK, , , 81 . A a 



b dK, , , , A „ 

M(*)*'-5-{!/}(«) =0 



dz 



dw 



for all x E [a, b], where K — XqL — AG. 

Corollary 16 (cf. Theorem 3.4 of |10j). Let y be a minimizer of 

J(y)= / L(x,y(x)CD«y(x))dx 
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subject to the boundary conditions 

y(a) = y a and y(b) = y b 

and the isoperimetric constraint 

,6 

I(y)= G{x,y{x), c a D«y{x))dx = 1 , 7 G R. 

J a 

Then, there exist two constants Xq and X, not both zero, such that y is a solution of equation 

^ (x, y(x), c a D«y{x)) + X D<£ (x, y(x), C a D«y{x))^ = 

for all x £ [a, b], where K = XqL — AG. Moreover, if y is not an extremal for I, then we may take 
A = l. 

5 Holonomic constraints 

In this section we consider the following problem. Minimize the functional 

J{y\,yi)= / L(x,yi(x),y 2 (x),^D^yi(x), a D^y 2 (x), a I^yi(x), a I^y 2 (x),z(x))dx, (23) 

J a 

where z is defined by 

z(x) = f X l(t, y 1 (t),y 2 (t)^D?y 1 (t)^Dfy 2 (t), a I t y 1 (t), a ^y 2 (t))dt, 

J a 

when restricted to the boundary conditions 

(yx{a),y 2 {a)) = and ( yi (b),y 2 (b)) = {y\,y b 2 ), yl,y a 2 ,y\,y\ 6 R, (24) 

and the holonomic constraint 

g(x,yx(x),y 2 (x)) = 0. (25) 

As usual, here 

(x,yi,y 2 ,v 1 ,v 2 ,w 1 ,w 2 , z) -> L(x, yi, y 2 , Vx, v 2 , Wx, w 2 , z), 
(x,yi,y 2 ,vi,v 2 ,wi,w 2 ) -> l(x,y 1 ,y 2 ,v 1 ,v 2 ,w 1 ,w 2 ) 

and 

(x,yi,y 2 ) -> g(x,yi,y 2 ) 
are all smooth. In what follows we make use of the operator [•, •] given by 

[j/i,2/2](a0 = (x, yi(x), y 2 (x),%D%yx(x), %D%y 2 (x), a lPyx(x), a I% y 2 (x), z{x)) , 

we denote (x, yi(x), y 2 (x)) by (x,y(x)), and the Euler-Lagrange equation obtained in ^ with 
respect to by (ELEi), i = 1,2. 

Remark 17. For simplicity, we are considering functionals depending only on two functions y\ 
and 2/2 • Theorem Q3 is, however, easily generalized for n variables yx, . . . , y n . 

Theorem 18. Let the pair (2/1,2/2) be a minimizer of J as in (|23[) , subject to the constraints 
(|24 [) — f|25|) . If 7^ 0, £/ien £/iere exists a continuous function A : [a, 6] — > R suc/i £/ia£ (2/1, 2/2) is a 
solution of 

(^fo>»K*)J+.tf (^'^V+I ^k.»K*)*^fa.»K») 

+ (^£^[yi,2/2](^-||-{2/i,2/2}(^+,/f (j"^[y 1 ,y 2 ](t)dt~ i {y u y 2 }(x)j=0 

(26) 

/or a; e [a, 6] and i = 1,2, where F[yx,y 2 ](x) = L[ yi ,y 2 ](x) - X(x)g(x,y(x)). 
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Proof. Consider a variation of the optimal solution of type 

(2/1,2/2) = (yi+£h 1 ,y 2 + eh 2 ), 

where h\,h 2 are two functions defined on [a, b] satisfying 

/11(a) =h 1 (b) = h 2 {a) = h 2 (b) = 0, 

and e is a sufficiently small real parameter. Since ■§^{x,y l {x),y 2 {x)) ^ for all x 6 [a, b], we can 
solve equation g{x,y 1 {x) 1 y 2 {x)) = with respect to h 2 , h 2 — h 2 (e,hi). Differentiating J(jJi,y 2 ) 
at e = 0, and proceeding similarly as done in the proof of Theorem [TJ we deduce that 



/ (ELE^h^x) + [ELE 2 )h 2 {x) dx = 0. 

J a 

Besides, since g(x,y 1 (x),y 2 (x)) = 0, differentiating at e = we get 



h 2 (x) 



Define the function A on [a, b] as 



X(x) 



&(*,y(*)) 
&(*,y(*)) 



(ELE 2 
Og 



hi{x). 



Combining (|28|) and ([29)) . equation (|27|) can be written as 



' \ELE 1 )-X(x)-^(x,y(x)) 



dyi 



h\(x) dx — 0. 



By the arbitrariness of hi, if follows that 



dg 



(ELE 1 )-X(x)-^(x,y{x))=0- 



Define F as 



F[y\,y 2 ]{x) = L[y l7 y 2 }(x) - X(x)g(x,y(x)). 
Then, equations (f2l)f follow. 



(27) 



(28) 



(29) 



□ 



6 Higher order Caputo derivatives 

In this section we consider fractional variational problems, when in presence of higher order Caputo 
derivatives. We will restricted ourselves to the case where the orders are non integer, since the 
integer case is already well studied in the literature (for a modern account see [T3lfT6l[29] ) . 

Let n £ N, (3 > 0, and € R be such that e (k — 1, k) for k e {1, . . . , n}. Admissible 
functions y belong to AC n ([a, b};R) and are such that ^D" k y, k = 1, . . . ,n, and y exist and 
are continuous on [a, b\. We denote such class of functions by T n {\a, b]; K). For a = (ati, . . . , a n ), 
define the vector 

c a D«y{x) = <°B?y(x), . . . °D^y{x)). (30) 
The optimization problem is the following: to minimize or maximize the functional 

J(y) = f L(x, y{x)°DZy{x), a I^y(x), z(x))dx, (31) 

J a 

y G J 771 ([a, bjjM), subject to the boundary conditions 

y {k) (a) = ya,k and y^(b)=y b>k , k € {0, . . . , n - 1}, (32) 
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where z : [a, b] — » R is defined by 

z{x) = f l(t,y(t)CD?y(t), a I?y(t))dt. 



Theorem 19. Ify G J- n (\a, &]; R) is a minimizer of J as in (|3ip . subject to the boundary conditions 
(|32p , fZien y is a solution of the fractional equation 

B i WW + g (s>m) + ^ (iww) + f ■ f <»> W 

+x>r(7 

k=l \ jx 



< /^ M (t)*.* M (,))-0 



/or a/Z a; £ [a, 6], wZiere = (x,y(x),^ D"y(x), a I^y(x), z(x)) with % D%y(x) as in ([50]). 

Proo/. Let Zi G J" n ([a, 6]; R) be such that ftW (a) = frW (6) = 0, for k G {0, . . . , n - 1}. Define the 
new function j as j(e) = J(y + eh). Then 



BT BT 8T 

g-[y](x)h(x) + ^2—[y}(x)^D^h(x) + —[y}{x)J^h{, 

k—1 



8L 
'ft?' 



/ x I ell " Pi] f)l \ 

l-{y}(t)h(t) + J2 ■j^{v}(t)ZB?'>h{t) + —{y}{t)J?h{t)\ dt 



dx = 0. (33) 



fntegrating by parts, we get that 



^\y]{x) c a D«*h[x)dx = x DZ*(J^[y}( X fjh(x)dx 



k-l 

E 



for all k G {1, . . . , n}. Moreover, one has 



b Q L 

Q^[y](x)al f xh(x)dx 



T ,3 ,'dL, 



J p b [-^[y]{x))h{x)dx, 



f "^[y](x)£^{y}(t)h(t)dtdx = J* n*^[y](t)dtj ^{y}(x 



)h(x) dx, 



dL 



dz [V][X) { L ^ k {y}i t ) C * D t k Ht)dt^ dx = J" ^ ^[y](t)d?j ^-mx) C a D^h[x) dx 



rb 



k-l 

E 



j-\aj,+m— k 



m=0 



b r)T f)l \ 

^[y](t)dt-^-{y}(x)\h^- m \x) 



-I b 



h QL dl \ 

[y](t)dt—{y}(x)}h(x)dx, 



dz 



dvk 
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and 



f a [[ ^{y}(t) a I?h(t)dt^ dx = jf\l£ ^ ^[y](t)dt^{y}{x)j h{x)dx. 

Replacing these last relations into equation p3[) . and applying the fundamental lemma of the 
calculus of variations, we obtain the intended necessary condition. □ 

We now consider the higher-order problem without the presence of boundary conditions Q32I) . 

Theorem 20. If y G J- n ([a, b]; K.) is a minimizer of J as in (|3ip . then y is a solution of the 
fractional equation 



dL 



E^r 



k=l 



dL 

dv k 



[v](x, 



dL 
dw^ 



i>L 'y](t)dt-^-{y}(x) 



dz 



dy 



for all x G [a,b], and satisfies the natural boundary conditions 



-i 6 



= 0, /or a// fee {!,..., n}. (34) 



Proof. The proof follows the same pattern as the proof of Theorem 1191 Since admissible functions 
y are not required to satisfy given boundary conditions, the variation functions h may take any 
value at the boundaries as well, and thus the condition 



h {k) (a) = h ik) (b) = 0, for k € {0, . . . , n - 1}, 



(35) 



is no longer imposed a priori. If we consider the first variation of J for variations h satisfying 
condition ()35|) , we obtain the Euler-Lagrange equation. Replacing it on the expression of the first 
variation, we conclude that 



fe-i 



EE 

k=l m=0 



a k +m-k 



0. 



To obtain the transversality condition with respect to fc, for fc £ {1, . . . , n}, we consider variations 
satisfying the condition 

fcCk-i)(o) ^ ^ ht k -V{b) and h^-^ia) = = h^~ l \b), for all j G {0, . . . ,n} \ {fc}. 



□ 

Remark 21. Some of the terms that appear in the natural boundary conditions (|34[) are equal to 
zero (cf. Remark^ and Remark XT^i . 



7 Fractional Lagrange problems 

We now prove a necessary optimality condition for a fractional Lagrange problem, when the 
Lagrangian depends again on an indefinite integral. Consider the cost functional defined by 

,6 

J{y,u) = / L (x,y(x), u(x) Igy(x),z(x))dx, (36) 
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to be minimized or maximized subject to the fractional dynamical system 

c a D«y{x) = f(x, y(x),u(x), a lPy(x), z{x)) (37) 
and the boundary conditions 

y(a) = y a and y(b) = y b , (38) 

where 



z(x)= / l(t,y(t)°D?y(t), a I?y(t) 



dt. 



We assume the functions (x,y, v, w, z) — )• f(x,y,v,w,z), (x,y,v,w, z) — > L(x,y,v,w, z), and 
(x, y, v, w) — > l(x, y, v, w), to be of class C 1 with respect to all their arguments. 

Remark 22. If f(x,y(x),u(x), a I@y(x), z(x)) = u(x), the Lagrange problem (|36|) -(|38 |) reduces to 
the fractional variational problem (j4|)-((5]) studied in Section^ 

An optimal solution is a pair of functions (y, u) that minimizes J as in (|36j) . subject to the 
fractional dynamic equation (|37j) and the boundary conditions fj38 j> . 



Theorem 23. If(y,u) is an optimal solution to the fractional Lagrange problem (136^ (|38[) . then 
there exists a function p for which the triplet (y,u,p) satisfies the Hamiltonian system 

OH 

aD°y(x) = —\y,u,p-](x), 

dU ( OH \ OH Qi 

x D%p{x) = — \y,u,p~}(x) + x l{} (— \y,u,p]{x)\ + J —\y,u,p](t)dt-—{y}(x) 

+xD » (rS" [y ' w ' pi(<) *'£ M(x) ) +xi > {i"^ iy ' u ' pmdt '§i {y}ix) 

and the stationary condition 

— \y,u,p\{x) = 0, 

where the Hamiltonian H is defined by 

H\y, u,p] (x) = L(x, y(x),u(x), a lgy(x), z{x)) + p(x)f(x, y(x),u(x), a lgy(x),z(x)) 

and 

\y,u,p](x) = (x,y(x),u(x), a lgy(x),z(x),p(x)), {y}{x) = (x,y(x),%D%y(x), a lPy(x)). 
Proof. The result follows applying Theorem [7] to 

,6 

J*(y,u,p)= / H\y,u,p](x)-p(x)%D%y(x)dx 

J a 

with respect to y, u and p. □ 

In the particular case when L does not depend on a I^y and z, we obtain [221 Theorem 3.5]. 
Corollary 24 (Theorem 3.5 of [22]). Let (y(x),u(x)) be a solution of 

r" 

J(y,u) = / L(x,y(x),u(x))dx — > min 



subject to the fractional control system ^D"y(x) = f(x,y(x),u(x)) and the boundary conditions 
y(a) = y a and y(b) = y b . Define the Hamiltonian by H(x,y,u,p) = L(x,y,u) + pf(x,y,u). Then 
there exists a function p for which the triplet (y, u,p) fulfill the Hamiltonian system 

{^D-y{x) = ^{x,y(x),u{x),p{x)), 
\ x D?p{x) = ^(x,y(x),u(x),p(x)), 

and the stationary condition ^-(x,y(x),u(x),p(x)) = 0. 
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8 Sufficient conditions of optimality 

To begin, let us recall the notions of convexity and concavity for C 1 functions of several variables. 

Definition 25. Given k £ {1, . . . ,n} and a function $ : D C R n — > K such that d^/dxi exist 
and are continuous for all i £ {k, . . . , n}, we say that "J is convex (concave) in (xk, ■ • ■ , x n ) if 



+ Ti, . . . ,X k -l + T k -l,X k + T k ,...,X n + T n ) - . . . ,X k -X,X k , ...,X n ) 



9* 

> (<) j. — {Xl, . . . , X k -l,X k , ■ ■ ■ ,X n )T k H 1- - — (xi, 

UXk OX n 



for all (xi, . . . , x n ), (x 1 +Tx,...,x n + r n ) £ D. 

Theorem 26. Consider the functional J as in and let y £ J-"([a, 6];R) be a solution of the 
fractional Euler- Lagrange equation ([6]) satisfying the boundary conditions ([5]). Assume that L is 
convex in (y, v, w, z). If one of the two following conditions is satisfied, 

1. I is convex in (y, v, w) and ^-[y](x) > for all x £ [a, b]; 

2. I is concave in (y, v, w) and ^\y\(x) < for all x £ [a, b]; 
then y is a (global) minimizer of problem (J4j — (J5j . 

Proof. Consider h of class F{[a, b];R) such that h(a) = h(b) = 0. Then, 
J(y + h)- J(y) 



> 



dL 



L[x,y(x) + h(x)^D^y(x)+^D-h(x), a I^y(x)+ a I^l 

l(t, y(t) + h(t)°D?y(t) + C a D«h{t), a I?y(t) + a I?h(t))<ttj dx 
L(x,y(x)CD«y(x), a lPy(x), f !(*,»(*),« D?y(t), jl ] y(t))dt)dx 



dx 



" V ^\y]{x)h{x) + ^\y]{x) c a D<Zh{x) + ^-[y](x) a I^h{x) 



dy 



dv 



dw 



^{y}(t)h(t) + ^- v {y}{t) c a D?h{t) + ^L{y}( t ) a I?h(t) ) dl 



dL 
dy 



MM 



(dL 
\ dv 



[*/](*) 



dL 
dw 



dL 



x 1 t 



f dL, ., , , dl r , 

Lir^-dw^ 



dl 



h{x)dx = 0. 

□ 



One can easily include the case when the boundary conditions ([5]) are not given. 



Theorem 27. Consider functional J as in (U) and let y £ J- ([a, b]} R) be a solution of the fractional 
Euler-Lagrange equation ^ and the fractional natural boundary condition (|17[) . Assume that L 
is convex in (y, v, w, z). If one of the two next conditions is satisfied, 

1. I is convex in [y, v, w) and ^[y](x) > for all x £ [a, 6]; 

2. I is concave in (y, v, w) and ^[y](x) < for all x £ [a, b]; 
then y is a (global) minimizer of ([?]). 
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9 Numerical simulations 



Solving a variational problem usually means solving Euler-Lagrange differential equations subject 
to some boundary conditions. It turns out that most fractional Eulcr-Lagrange equations cannot 
be solved analytically. Therefore, in practical terms, numerical methods need to be developed and 
used in order to solve the fractional variational problems. A numerical scheme to solve fractional 
Lagrange problems has been presented in [1] . The method is based on approximating the problem 
to a set of algebraic equations using some basis functions. A more general approach can be 
found in [Uj] that uses the Oustaloup recursive approximation of the fractional derivative, and 
reduces the problem to an integer order (classical) optimal control problem. A similar approach is 
presented in 25 , using an expansion formula for the left Riemann-Liouville fractional derivative 
developed in [TT] , Here we use a modified approximation (see Remark [29f based on the same 
expansion, to reduce a given fractional problem to a classical one. The expansion formula is given 
in the following lemma. 

Lemma 28 (cf. equation (12) of pi]). Suppose that f E AC 2 [0,b], f" € Li[0,6] and < a < 1. 
Then the left Riemann-Liouville fractional derivative can be expanded as 

oo 

D%f(x) = A(a)x- a f{x) + Bi^x^f'ix) - ^ C(k, ^^-"^(i), 



k=2 



where 



v'k{x) 
A(a) 

B(a) 
C(k,a) 



= (l-k)x k - 2 f(x), v k (0)=0, fc = 2,3,.. 

T(k -l + a) 



r(l-a) r(2 - a)T(a - 1) E (fc-1)! ' 



k=2 



r(2-a) 



k=l 



1 



T(k - 1 + a) 
T(a- l)fc! 

T(k- 1 + a) 



r(2 - a)T(a - 1) {k - 1)! 
In practice, we only keep a finite number of terms in the series. We use the approximation 

N 



D2f(x) ~ A{a, N)x~ a f{x) + B{a, N)x l - a f{x) - ^ C(k, a)x 1 - k - a v k (x) (39) 
for some fixed number N, where 



fc=2 



N 



A(a,N) 
B{a,N) 



r(l-a) r(2 - a)T(a - 1) ^ (fc-1)! 



r(fc - l + a) 



1 



r(2-a) 



JV 



fc=i 



r(fc - 1 + a) 
T(a- l)fc! 



Remark 29. In /jij / i/ie authors use the fact that 1 + X^fcli ^(L-i)^ = ®> an< ^ a PP^V * n ^ le * r 
method the approximation 



N 



D a x f{x) ~ i4(a, N)x~ a f{x) - £ C(fc, a)* 1 -*""^^ 



fc=2 



Regarding the value of B{a,N) for some values of N (see Ta&ZeQJ), we take the first derivative in 
the expansion into account and keep the approximation in the form of equation (|39|) . 
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N 


4 


7 


15 


30 


70 


120 


170 


B(0.3, TV) 


0.1357 


0.0928 


0.0549 


0.0339 


0.0188 


0.0129 


0.0101 


B(0.5, N) 


0.3085 


0.2364 


0.1630 


0.1157 


0.0760 


0.0581 


0.0488 


B(0.7, N) 


0.5519 


0.4717 


0.3783 


0.3083 


0.2396 


0.2040 


0.1838 


5(0.9, N) 


0.8470 


0.8046 


0.7481 


0.6990 


0.6428 


0.6092 


0.5884 



Table 1: Values of B(a,N) for a e {0.3,0.5,0.7,0.9} and different values of N. 



We illustrate with Examples [3] and 0] how the approximation (131)1) provides an accurate and 
efficient numerical method to solve fractional variational problems. 

Example 30. We obtain an approximated solution to the problem considered in Example^ Since 
y(0) — 0, the Caputo derivative coincides with the Riemann-Liouville derivative and we can 
approximate the fractional problem using (|39p . We reformulate the problem using the Hamiltonian 
formalism by letting g 7 D"y(x) — u(x). Then, 



N 



A(a, N)x~ a y{x) + B(a, N)x 1 - a y'{x) - a Y- 



,1-fc- 



"v k (x) = u{x). 



(40) 



k=1 



We also include the variable z(x) with 

z'{x)= (y{x)-x a+1 )\ 
In summary, one has the following Lagrange problem: 

J{y)= [ [(u(x) -T{a + 2)x) 2 + z(x)]dx 



'y'(x) = -AB^x^yix) + ££L 2 B^Cux^v^x) + B- l x a - l u{x) 
v' k {x) = (l-k)x k - 2 y{x), k=l,2,... 
y(x) = {y(x)-x<*+ 1 ) 2 



(41) 



subject to the boundary conditions y(0) = 0, z(0) = and Ufc(0) = 0, k = 1,2,... Setting N = 2, 
the Hamiltonian is given by 

H = -[(u(x) -T(a + 2)x) 2 + z{x)]+p x {x) {-AB- 1 x~ 1 y{x) + B^C^v^x) + B^x^uix)) 

-p 2 (x)y(x) +p 3 (x) (y(x) - x a+1 ) 2 . 

Using the classical necessary optimality condition for problem (]4ip . we end up with the following 
two point boundary value problem: 



f y'W 



-AB- X x- X y{x) + B- 1 C 2 x- 2 v 2 (x) + \B- 2 x 2a - 2 p x {x) + T(a + 2)B- 1 x a 



v' 2 {x) = -y{x) 

z'(x) = (y(x) -x a+1 ) 2 

p[(x) = AB- l x- l px{x)+p 2 {x) -2p 3 (x)(y(x) - x a+1 ) 



(42) 



P 2 (x) 



-B^C^piix) 



1 



subject to the boundary conditions 



y(0) = 
«a(0) = 
1 2(0) = 



1/(1) = 1 
Pa(l) = 

b 3 (i) = o. 



(43) 



We solved system (|42|) subject to (I43|) using the MATLAB® built-in function bvp^c. The resulting 
graph for y(x), together with the corresponding value of J , is given in Figured 
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0.2 0.4 0.6 0.8 1 



x 

Figure 1: Analytic vs numerical solution to problem of Example [3] 



Our numerical method works well, even in the case the minimizer is not a Lipschitz function. 

Example 31. An approximated solution to the problem considered in Example^ can be obtained 
following exactly the same steps as in Example \3(A Recall that the minimizer to that problem 
is not a Lipschitz function. As before, one has y(0) = and the Caputo derivative coincides 
with the Riemann-Liouville derivative. We approximate the fractional problem using (|39j) . Let 
qD"u(x) = u(x). Then (|40[) holds. In this case the variable z(x) satisfies 

z'(x) = ( y(x) - ^ 

w V r(a + i)y 

and we approximate the fractional variational problem with the following classical one: 

J{v) = I [( u ( x ) ~ I) 2 + z(x)]dx — > min 
Jo 

'y'{x) = -AB- 1 x- 1 y(x) + J2k=2 B ' 1C k x ' k Vk(x)+B- 1 x a - 1 u{x) 
v' k (x) = (l-k)x k - 2 y(x), fc = l,2,... 

z'{x) = (y(x) - rgpy) 



2 



subject to the boundary conditions y(0) = 0, z(0) = and Ufe(0) = 0, k = 1, 2, . . . Setting N = 2, 
the Hamiltonian is given by 

H = -[(u(x) - l) 2 + z(x)] +pi(x) (—AB~ 1 x~ 1 y(x) + B' 1 C 2 x- 2 v 2 {x) + B~ 1 x a ~ 1 u(x)) 

-p 2 (x)y(x) +p 3 (x) (y(x) - T (* + 1 ^ • 
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The classical theory \36^ tell us to solve the system 

' y'(x) = -AB- 1 x- 1 y{x) + B~ 1 C 2 x~ 2 v 2 (x) + \B- 2 x 2a - 2 p x (x) + B~ 1 x c 

v' 2 {x) = -y(x) 

z\x) = (y(x) - p^) 2 

p[(x) = AB- l x-^{x) + p 2 (x) - 2p 3 (x)(y(x) - j^f^y) 

p' 2 (x) = -B- 1 C 2 x- 2 pi(x) 



(44) 



subject to boundary conditions 



'y(0) = o 

wa(0) = 
z(Q) = 



P2(l)=0 

,p 3 (l) = 0. 



(45) 



As done in Examvle \3(A we solved (I44[) - (|45[) using the MATLAB® built-in function bvp4c. The 
resulting graph for y{x) , together with the corresponding value of J, is given in Figure\^in contrast 
with the exact minimizer 



1.4 



1.2 



■ Analytic: J=0 

Approximation: N=2, J=0. 00042332 




0.4 0.6 

x 



Figure 2: Analytic vs numerical solution to problem of Example |U 
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